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Exact microscopic spectral correlation functions are derived by means of the replica limit of the 
Toda lattice equation. We consider both Hermitian and non-Hermitian theories in the Wigner- 
Dyson universality class (class A) and in the chiral universality class (class AIII). In the Hermitian 
case we rederive two-point correlation functions for class A and class AIII as well as several one- 
point correlation functions in class AIII. In the non-Hermitian case the average spectral density of 
non-Hermitian complex random matrices in the weak non-Hermiticity limit is obtained directly from 
the replica limit of the Toda lattice equation. In the case of class A, this result describes the spectral 
density of a disordered system in a constant imaginary vector potential (the Hatano-Nelson model) 
which is known from earlier work. New results are obtained for the average spectral density in the 
weak non-Hermiticity limit of a quenched chiral random matrix model at nonzero chemical potential. 
These results apply to the ergodic or e domain of the quenched QCD partition function at nonzero 
chemical potential. Our results have been checked against numerical results obtained from a large 
ensemble of random matrices. The spectral density obtained is different from the result derived 
by Akemann for a closely related model, which is given by the leading order asymptotic expansion 
of our result. In all cases, the replica limit of the Toda lattice equation explains the factorization 
of spectral one- and two-point functions into a product of a bosonic (noncompact integral) and a 
fermionic (compact integral) partition function. We conclude that the fermionic partition functions, 
the bosonic partition functions and the supersymmetric partition function are all part of a single 
integrable hierarchy. This is the reason that it is possible to obtain the supersymmetric partition 
function, and its derivatives, from the replica limit of the Toda lattice equation. 



I. INTRODUCTION 



A replicated system is one in which certain fields appear as exact copies of each other. For a given number, say n, of 
replica fields the system may or may not have a physical realization, but it is nevertheless useful to study the system 
for any integer values of n. Suppose, as is the case in this paper, that we wish to calculate the average resolvent of 
some operator D, 



and, as is also the case, that it is simpler to evaluate (det n (D + x)) for integer values of n than it is to evaluate the 
average of the resolvent. In that case the replica trick allows us to obtain the resolvent from the replicated system. 
The fermionic (bosonic) replica trick consists in evaluating (det™(D + x)) for all positive (negative) integers, n, and 
then taking the limit n^Oas follows [1] 



The tricky part arises because the use of the replica trick assumes that an analytic continuation from positive or 
negative integers is uniquely defined. It has long been known [2-12] that this is the case for perturbative calculations 
such as for the expansion of the resolvent in an asymptotic series in 1/x. Problems are known to occur when one 
attempts to derive exact nonperturbative results for the quenched average resolvent [4,7]. This is one of the reasons 
that alternative methods were introduced to evaluate quenched averages. Most notably, the supersymmetric method 
[15,16], where the resolvent is obtained from the generating function 



has been very successful. Among other alternatives to the replica trick, we only mention the Keldysh method [17] by 
means of which some nonperturbative results have been derived. 

The purpose of this paper is to show that nonperturbative results can be obtained by means of the replica trick as 
well. 

A determinant raised to an integer power inside a statistical average appears in a wide range of theories. In 
disordered condensed matter systems and random matrix theory, one is usually interested in quenched disorder where 
such a determinant enters in the generating function of Greens functions as discussed in our example of the resolvent. 
In Quantum Chromodynamics Dynamics (QCD), on the other hand, the physical theory is an average over gauge 
field configurations of a fermionic determinant raised to the number of flavors. In lattice gauge simulations, because 
of computational reasons, one often is forced to work with quenched averages (i.e. ignoring the fermion determinant) 
which arc not realized in nature. Both in disordered systems and in QCD the low-energy effective theory of the 
quenched system is either a supersymmetric nonlinear sigma model or a nonlinear sigma model for replicated flavors. 
For perturbative calculations, both theories are completely equivalent [4,8,11-14]. To prefer one above the other is 
just a matter of convenience. On the other hand, nonperturbative results have only been derived by means of the 
supersymmetric method [15,16,18-22]. 

The replica trick has met substantial critique [4,7] in nonperturbative applications. For instance, it was shown 
that bosonic and fermionic replica limits lead to different results for the two point function of the Gaussian unitary 
ensemble and that neither limits lead to the correct result [4]. The original critique [4] motivated an insightful and 
challenging debate [5-10,23-25] the details of which have been accounted for in the introduction of [26]. The debate 
lead to the generally accepted belief that, although the replica trick correctly reproduces the asymptotic expansion of 
the microscopic spectral correlation functions, it can not be trusted for nonperturbative calculations except in a few 
cases, when the series can be resummed [24] or simply terminates [7,9]. However, generally this is not the case, and 
the asymptotic series leaves much to be desired. For example, since an asymptotic series does not determine the cut 
structure of an analytic function, it is not possible to obtain the density of eigenvalues from the discontinuity of the 
resolvent across its cut. 

Over the past few years intimate connections have been established between random matrix theory, integrable 
systems, the theory of r-functions, Painleve equations and unitary matrix integrals [27-31,35,8,36,32-34]. Kanzieper 
[37] exploited this connection to obtain exact analytical results by means of the replica limit of Painleve equations 
for the generating function of spectral correlation functions. He noted that the replica index only appears as a 
coefficient in the associated Painleve equation. To leading order in the replica index this equation is solvable and, 
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with appropriate boundary conditions, exact analytical results for spectral one- and two-point functions functions 
could be derived. 

The ground breaking work of Kanzieper was followed up by the present authors [26] by showing that the replica 
limit can be taken in the Toda lattice equation rather than in the Painleve equation. This method has two advantages. 
First, it automatically reveals the factorization of correlation functions or derivatives of correlation functions in terms 
of a product of a fcrmionic and a bosonic partition function. Second, in the examples we have considered, cither 
no differential equation had to be solved or the desired result followed from the solution of a simple first order 
(inhomogeneous) differential equation. The factorization property, which had not been appreciated before, explains 
precisely why the original formulation of the replica trick is bound to be tricky: Fcrmionic partition functions are 
given by compact integrals which are analytic functions of the argument for all integer values of the replica index 
n. It is therefore a mystery how the replica limit, n — > 0, can be a nonanalytic function. This point of the original 
critique [7] is dealt with automatically in the replica limit of the Toda lattice equation which involves both fcrmionic 
and bosonic partition functions. The bosonic partition function is determined by a noncompact integral and carries 
the nonanalytic information. 

The Toda lattice is well-known from the theory of exactly solvable systems [38,39]. It consists of masses interacting 
in one dimension through an exponential potential [38]. The Hamilton equations of motion are known as the Toda 
lattice equation and its solutions give the position of the masses as a function of time. The system is intcgrable 
and one particular solution can be written as a unitary matrix integral which, among others, is the QCD partition 
function in the e regime. Similar connections exist for all of the classical random matrix ensembles (see for example 
[40]). Remarkably, the replica trick also enters in integrable systems related to matrix models for supersymmetric 
Yang-Mills theories [41], but it is not known if there is a relation with our approach. 

In [26] we derived the two-point function of the Gaussian unitary ensemble and the quenched and partially quenched 
one-point functions of the chiral Gaussian unitary ensemble by means of the replica limit of the Toda lattice equation. 
The purpose of the present paper is to show that such strikingly simple results apply more generally. We will consider 
three different physical systems: the QCD partition function, the Hatano-Nelson model [42] and QCD at nonzero 
chemical potential. The Hatano-Nelson model is a disordered system in an imaginary vector potential that enters in the 
same way as a baryon chemical potential in the Euclidean QCD partition function. In both systems, the symmetries 
of the generating function are broken spontaneously. The corresponding Goldstone bosons interact according to a 
nonlinear sigma model that is completely determined by the pattern of spontaneous and explicit symmetry breaking 
[43]. In QCD, this is the well-know chiral Lagrangian which can be extended to include a nonzero chemical potential 
[44]. For the Hatano-Nelson model a similar effective Lagrangian can be derived [15]. Its static limit was obtained 
earlier [45] as the weak non-Hermiticity limit of a complex non-Hcrmitian random matrix model. 

In both theories we will study the nonlinear sigma model in the ergodic domain which, in QCD, is also known as the 
e-regime, where the contribution from the zero-momentum modes (or constant fields) dominate the partition function. 
In this domain, the kinetic term of the a model can be ignored. The corresponding energy scale, which is known as 
the Thouless energy, is the quark mass for which the Compton wave length of the Goldstone bosons is of the order 
of the size of the box. In this regime, all theories with a given symmetry breaking pattern are universally described 
by the same partition function. Perhaps the simplest theories in this universality class are random matrix theories. 
In random matrix theories one seeks to describe average spectral properties on the scale of the average eigenvalue 
spacing in the limit of large matrices. It has been shown that spectral correlation functions in this microscopic limit 
are universal in the sense that they only depend on the symmetries of the random matrix model and not on its 
detailed form (see for example [46,47]). We now understand that [21,20] this universality is a direct consequence of 
the uniqueness of the static part of the effective Lagrangian. We will exploit the correspondence between random 
matrix theory and the static part of the effective Lagrangian to derive the bosonic partition partition function for 
QCD at nonzero chemical potential and for the Hatano-Nelson model. In both cases, the integration contour of the 
bosonic fields in the cr-model follows naturally from the Ingham-Siegal integral of the second kind [25]. 

The two-point functions of the unitary ensemble, the chiral unitary ensemble and the spectral density of the 
Hatano-Nelson model [45,15] and the chiral unitary ensemble at nonzero chemical potential all share one remarkable 
property. They factorize into a product of a fcrmionic and a bosonic partition function. We will show that this is a 
natural consequence of the Toda lattice equation. Indeed, the corresponding generating functions can be written as 
a determinant of derivatives. Therefore, they are r-functions of a Toda lattice hierarchy and satisfy a Toda lattice 
equation. Our result for the chiral unitary ensemble at nonzero chemical potential is a genuinely new result which 
is in complete agreement with results obtained from a numerical matrix diagonalization. The leading term of an 
asymptotic expansion of our result and the result for a closely related model for non-Hermitian spectra with a chiral 
symmetry [48] turn out to be the same. 

The organization of this paper is as follows. As a warm up we start in sections II A and II B with simple examples 
that have already been discussed in [26]. The two-point function of the chiral unitary ensemble is discussed in 
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section II C. Turning to the complex spectra of non-Hcrmitian operators we consider the case of a disordered system 
in a nonzero imaginary vector potential in section III A and QCD at non zero chemical potential in section IIIB. 
Concluding remarks are made in section IV. Our notation as well as additional technical details are explained in 
several appendices (A through D). 



II. TODA LATTICE EQUATION FOR HERMITIAN THEORIES 

In this section we start with three simple examples of the replica limit of the Toda lattice equation, i.e. the two- 
point function of the Unitary Ensemble (UE), the resolvent of the chiral Unitary Ensemble (chUE) in the quenched 
case for topological charge v, and in the case of Nf massless flavors in the sector of zero topological charge. The first 
two cases were already discussed in [26] . In the second part of this section we discuss the calculation of the two-point 
function of the chUE. In this case the Toda lattice equation is obtained from a new representation of the finite volume 
partition function which might also be useful in other applications. 



A. The Two-Point Function of the Unitary Ensemble 

This section discusses the simplest example of the replica limit of the Toda lattice equation [26], namely the 
connected two-point function of the Gaussian Unitary Ensemble (GUE). Its generating function for n fermionic 
replicas is given by 

J dH dct n (H + E±)dct n (H + E 2 )e-% TtH>H . (4) 

Following the usual reduction, this partition function can be written as a nonlinear cr-model. In the large iV-limit, at 
fixed r — N(E2 — E\) near the center of the spectrum, the eigenvalue representation of this cr-model further simplifies 
to (see eq. (1.2) of [4]) 

Z n (ir) = n\ f f[ d\ k e-* rX *A 2 (\). (5) 
fe=i 



Here, A(A) is the Vandermonde determinant A(A) = Yl k<t (^k — A;). The connected two-point correlation function is 
then obtained simply as the replica limit 



G{r) = - lim \d 2 r log Z n (ir). (6) 



The generating function (5) can be rewritten as 

Z n (r) = (n!) 2 det[a* + ^i(r)] w=0> ...,„-i. (7) 
In a similar way the large TV limit of the bosonic partition can be reduced to 



i r 00 _ n _ 

Z -^) = ^ W J i n^A* ( A). 



(8) 



where the sign of the exponent is determined by the sign of the imaginary part of r. The r dependence can be scaled 
outside the integral resulting in 



Z -n(r) = ^U(^ 2 - ( 9 ) 
fe=i 

This expression can be rewritten as 

Z - n{r) = [(n-l)!]* det $ +Ji Mr)k;=o,..,n-i- (10) 
Therefore, both the fermionic and the bosonic partition function satisfy the Toda lattice equation [23,26] 
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- z U ir )dr log Z n {ir) = j-—^Z n+1 {ir)Z n - 1 {ir). (11) 

In the replica limit, the l.h.s. of this equation is n 2 times the connected two-point function whereas the r.h.s. correctly 
gives 

G(r) = Z x {ir)Z_ x {ir) = ^EI!^. (12) 

We thus find that the factorization of the two-point correlation function into a bosonic and a fermionic partition 
function is not accidental but rather a consequence of the relation between random matrix theories and integrablc 
hierarchies. The two-point function of the GUE can also be obtained from a solution of the Painleve V equation [23]. 
Another feature of this result emphasized in [7] is that the asymptotic series in 1/r terminates. This explains [7] that 
in this case the exact result could be obtained [5] by including a sublcading saddle point manifold in the standard 
replica approach. 



B. The Partially Quenched Resolvent in QCD 

In this section we discuss the replica limit of the Toda lattice equation for the generating function of the partially 
quenched resolvent in chUE. As a simple extension of the results presented in [26] we obtain the microscopic limit of 
the partially quenched resolvent for an arbitrary number of flavors. 

The partition function of the chiral Gaussian Unitary Ensemble (chGUE) with n fermionic replica flavors with mass 
m and Nf additional fermionic flavors with masses mi, ■ ■ ■ , mjv, is defined by 

N f 

Z n,N f (m, {mi}) = (det n (D + m) JJ det(L> + m k )). (13) 

k=i 

The Dirac matrix is given by 

iW 



D =liWt )' ^ 



where W is a rectangular lx(l + v) matrix so that D has exactly v zero eigenvalues. The average is over the probability 
distribution of the matrix elements of W with the weight 

P(W) = e -i N ^ wwt \ (15) 

with N = 21 + v. We will consider the microscopic limit where the thermodynamic limit is taken for fixed x = mJV 
and Xk = rrikN. In this case the partition function reduces to the unitary matrix integral [49] 

Z {v) (xtx-\) = - f dUdet"Uei Tr[Miu+Mui] (16) 

where M — diag(ir, • • • , x, x\, ■ ■ ■ , Xn.) is the mass matrix. The normalization factor has been included for later 
convenience. The partially quenched resolvent, which we will derive below, is given by 

G(x,{ Xi }) - hm i -d x logZi% Ax,{ Xi }). (17) 

n— >0 n ' 1 

The unitary matrix integral (16) is the partition function of QCD [49-51] in the e-regime where it coincides with the 
microscopic limit of the chGUE partition function (13). Partition functions for different number of flavors satisfy a 
recursive relation know as the Toda lattice equation [30] 

E f a a i ryM , r -, s 2 Z n+1, Nf ( X > { X k }) Z n- 1, N f ( X > { X k } ) 
xod Xo x k d Xk logZ"> N (x ,{xk}) =nx — — — . (18) 

For degenerate masses x k the sum over the derivatives is only over different masses. 
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The bosonic chGUE partition function is defined as in (13) but with negative integer values of n. Its axial symmetry 
is not U(n + N f ) but rather Gl(n + N f )/U(n + N f ) [20,21,9]. Therefore, the Goldstone manifold for the bosonic 
partially quenched partition function is Gl(n + Nf)/U(n + Nf) rather than U(n + Nf), i.e. the coset of positive 
definite matrices. In this way the integration manifold remains Riemannian [20] so that all integrals are convergent 
for positive masses. As bosonic partition function we thus obtain [21,9] 



If we parameterize Q as 
the invariant measure is given by 



(x,{xi}) = ^-[ dQde^Qe-^^Q+^l (19) 

k-n JQeGl(n)/U(n) 

Q = Udia,g(e Sk )U\ U^U = 1 (20) 



dQ = ]J(e Sk -e s ')(e- Sk - e~ si ) J| ds k dU. (21) 



k<l 



Using the Ingham-Siegel integral [25] this result can be derived directly from the bosonic random matrix partition 
function [52]. For zero topological charge this is discussed in [25] and the simplest nontrivial example is given in 
Appendix C. The explicit form of the partition function with any number of fcrmions and bosons in the topological 
sector v conjectured in [26] was proven rigorously in [53] based on a general result in [54]. 

In [26], the n — > limit of the Toda lattice equation (18) was solved for Nf = and Nf = 1. Before we consider 
the case with Nf massless flavors and zero topological charge, we first review the quenched limit (Nf = 0) of (18) in 
the sector of topological charge v. In this case, (18) simplifies to 

[xd x f log Z W (x) = nx 2 {X) - (22) 

[Z n (x)\ 

Collecting terms of 0(n) in the replica limit of this Toda lattice equation we find 

d x [xG(x)} =xZ^\x)Z ( l\(x). (23) 

Again, we observe a factorization property of the resolvent. For topological charge v we have that * 

Z { "\x)=I v (x) and Z { l\(x)=2K v (x). (24) 

The general solution of the inhomogeneous differential equation (23) is given by 

xG(x) = a + x 2 (K u (x)I v (x) + K^^I^x)), (25) 

where the first term is a solution of the homogeneous differential equation. The integration constant is fixed by the 
small x expansion of the resolvent which in the sector of topological charge v is given by 

G(x) ~ -, (26) 
x 

so that a = v. This ends the derivation of the quenched resolvent in the topological sector v. The result is in exact 
agreement with the known result [55]. 

As a final introductory example, let us now consider the case of Nf massless flavors and zero topological charge. 
The general Toda lattice equation (18) reduces to 



1 



lim ^ (xd x + yd y )xd x log Z { ^ Nf (x, y) 



= x 2 Z ^-,": (27) 

[z { ^ s (xM 2 



%(x,0)Z% ; (x,0) 



*The normalization factor of the bosonic partition function is chosen such that the large x limit of both sides of (23) is the 
same. For large x we have that G(x) — 1 whereas the large x behavior of the r.h.s. of (23) is obtained from the leading order 
asymptotic expansion of the Bessel functions. 
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The partition function Z\ ' N (x, 0) can be obtained by using the flavor-topology duality relation [56] and is given by 

Z[% f (x,0) = x- N 'Z[j N '\x) = x- N 'I Nf {x). (28) 

The partition function Z^\ N (x, 0) has Nf massless fermionic flavors and one bosonic flavor with mass x. Using a 
bosonic version of the flavor-topology duality relation we find 

Z { Xn s {x, 0) = x N fZ% Nf \x) = 2x N fK N; (x), (29) 

where the normalization factors are consistent with the asymptotic behavior of (27). The y<9 y -derivative of 
log Z^Nf i x ' V) vanishes for y = so that (18) simplifies to 

8 X [xG{x)\ = 2xK Nf (x)I Nf (x) . (30) 

As in the quenched case for topological charge v, the solution is given by 

xG(x) = a + x 2 [K Nf (x)I Nf (x) + K Nf ^ 1 (x)I Nf+1 (x)}, 

Because we are in the sector of zero topological charge, the resolvent approaches a constant for small x so that the 
integration constant a = 0. This result agrees with the result obtained by integration the microscopic spectral density 
[55]. 



C. The Two-Point Function in QCD 

In this section we consider the microscopic limit of the chUE partition function for m flavors with mass x and n flavors 
with mass y. In the sector of topological charge v we will denote this partition function by Zm]n{x, y). The microscopic 
spectral two-point function will be derived from the replica limit of the Toda lattice equation corresponding to this 
partition function. In the first subsection we will show that the partition function can be written as a determinant 

of derivatives. In the second subsection we shall use this form to show that Zm]n{x,y) satisfies a Toda lattice 
equation and therefore is a r-function of an integrable hierarchy. The replica limit of the Toda lattice equation then 
automatically gives us the exact analytical result for the two-point function that was obtained earlier [57] by means 
of the supersymmetric method. 



1. The Generating Function for the Two-Point Function is a r— Function 

The finite volume partition function for m fermionic flavors with mass x and n (with n>m) fermionic flavors with 
mass y is given by [49-51] 

This partition function can be derived from a random matrix theory with the symmetries of the QCD partition function 
[49] and is also known as the QCD partition function in the e regime. This integral and its generalizations have been 

well studied in the literature [58-62], but the known analytical expressions do not directly show that Zm]n{x,y) is 
a r-function. In order to write it in such a form we evaluate the integral over U by decomposing U(n + m) as 
U(m) x U(n) x [U(n + m)/(U(n) x U(m))], and evaluating the integral over the coset first. The normalization factor 
is chosen to be 

T m , n = (2 7r )(«+ m )(™+ m + 1 )/ 2 . (32) 
Our calculations simplify if we use the parameterization 
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where A is the block diagonal (m + n) x (m + n) matrix given by 



A = ( /i -v^~^ | • (34) 
0-1 



The diagonal matrix fi — diag(/xi, • • • , /i m ), and we will use the notation that A^ = yl— 7*1 with A& G [0, 1]. The 
unit matrix 1 in the lower right hand corner is of size n — m. Furthermore, u\ and v\ are unitary m x m matrices, 
i*2 € U(n) and v 2 ^ U (n)/(U m (l) x U(n — m)). One easily verifies that the total number of parameters on both sides 
of (33) is the same. The Jacobian of this transformation is derived in Appendix B and is given by 

m 

J= II (A 2 fe -Af) 2 n(2A*)M 2 fe ( "- m) . (35) 

l<k<l<m fe=l 

The finite volume partition function can thus be rewritten as (the Vk variables have been eliminated from the exponent 
by the transformation u k — * VkU^v^ 1 ) 

z $n( x >v) = — tL / dv i [ dv 2 [ dm [ du 2 ( f\d\ k 

m\± m ,n JU(m) J U(n) /[U m (l)xU(n-m)] JU(m) JU(n) JO k=1 

x J({Afe})dct ,y (u lU2 )exp[ixTrA 11 (u 1 + u\) + ^TrA 22 (u 2 + 4)]. (36) 
Here we have introduced the notation 

An = diag(Ai, . . . , A m ), (37) 

and 

A 22 = diag(-Ai, . . . , -A m , -1, . . . , -1). (38) 

The integrals over V\ and w 2 give an overall constant equal to the volume of the integration manifold. The integrals 
over Mi and w 2 are well known. They are given by [59-62] 



1 f 7 i \ rrr-i / A , 1 A fM n d0t ^(^fc) fc,i=0,..., m -l 

77777 ^ / a^idet (m) exp Tr(-a;Aiitti + -iAn« = C m . 

vol [/ m ))J u(m) 2 2 A { xA fe 2 } 



(39) 



where 

<t>l{x k ) = (x k d Xk ) l I v {x k ), (40) 

and 

m— 1 

C m = 2 m < m - 1 )/ 2 J| fc! . (41) 
fc=i 

The volume of the unitary group is given by (see eg. [63]) 

( 27r )n(n+l)/2 

W) life I 1 fc! 



vol(C/(n)) = / d£/ = V ■ (42) 



A similar formula can be derived for the u 2 -integral. The only complication is that n — m diagonal matrix elements 
of j/A 22 are degenerate. To apply the formula (39) we add e k to the degenerate diagonal matrix elements. The e k — > 
limit can then be obtained conveniently by writing the Taylor expansion of <pi{—y + e k ) to order n — m in e k as a 
product of a matrix containing the powers of e k and a derivative matrix. This results in 
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1 f 11 

UTU \ \ / du 2 det 1/ (u 2 )exp[Tr(-2/A22U2 + -yk 22 u\)] 
vol(U(n)) J V ( n ) Z l 



fi in— m)(n— m— l)/2 

_^ I det 

Cn-m TYL^IV 2 V 2 ) {n - m) ni< fe<i < m (Aly 2 - \W) 



<j>o{-yX\) ••• ^„_i(-j/Ai) 

0o(-yA m ) ••• <j> n -i{-yX m ) 
<f>o(-y) ■■■ <p n -i(-y) 



^- m_1) (-») ••• fc m_1) (-») 



(43) 



where we have used the notation (f>^\y) = d y 4>k(y)- All factors in the Jacobian (35) except IK^A/c) are canceled by 
the corresponding factors from the u\ and u 2 integrations. Multiplying the two determinants the partition function 
can be written as 



t){m+n)(m+n — l)/2„.(n— m)(n — m— l)/2 

7( v ) ( T v) — - - det 

^m,n^.W CnCmX rn( m -l) y n(n-l) 



d\\(j) (\x)(j) (-\y) ■■■ Jq 1 dAA0o(Aa;)0„-i(-Ay) 
Jq 1 d\\(j) m -i{\x)^Q{-\y) ■■■ ^ d\\(j) m -i{\x)(t) n -i{-\y) 



<t>o(-y) 



An-m-l). 



4>n-i(-y) 



i n -r 1} (-y) 



By the addition of rows, the derivatives in the last m rows can be rewritten in terms of the derivative operators 

d 



• (44) 



S x — x 



dx 



(45) 



This results in the final form for the partition function of the chUE 

£ d\\I v {\x)I v {-\y) 



z% n {x,y) 



2(m+n)(m+n— 1)/2 
C n C m X m » ~ 1 )y n ( n ~ !) 



det 



S^ 1 J^dXXI v (Xx)I v (-Xy) 



8™- 1 ft dXXI v {Xx)I v {-Xy) ■■■ ajT 1 ^" 1 £ dXXI v {Xx)I v {-Xy) 



lxn-l f 1 . 



iA-y) 
s^-'iA-y) 



5l n - m - 2 i v {- y ) 



(46) 



If we denote the n x n matrix that enters in this partition function by A^' m ^ it can be rewritten as 

a(m+n)(m+n— 1)/2 

Z^ n (x,y) = — — , . detA^ n \x,y). 

' C n C m x m ( m -Vy n( - n ~^ 



(47) 



Of course, A^ m ' n>> is a function of v even though we have not explicitly indicated it. Before we proceed to show that 
this partition function satisfies a Toda lattice equation, let us note that: z[^(x,y) —2j dXXI v {Xx)I v {—Xy). Hence, 
for tl — tyl we obtain the compact expression 



Another special case is m — 0, where the partition function simplifies to 

i»(-y) ••• S^I^-y) 
— — — det : 

iy v 



(48) 



2 n(n-l)/2 
Z m-0 n( X >y) = 7 TT ^t 



(49) 



5^I v {-y) ■■■ S 2 y n - 2 L{-y) 
This expression can be easily rewritten in terms of one of the familiar forms for the finite volume partition function. 
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2. The Toda Lattice Equation 

From the structure of Zm]n{x,tj) in (46) it automatically follows that Zm]n{x,y) satisfies a Toda lattice equation. 
Below, this will be shown by a slight extension of an argument given in [40] . 

For the determinant of a matrix A we have the Sylvester idendity [64] 

CijCpq CiqCpj = det (A) Cij t pq (50) 

where Cij is the cofactor of matrix element ij 

and Cij^p q is the double cofactor of matrix elements ij and pq 

_ d 2 det(A) 

C ™i ~ dA^dApq- (b2) 

We apply this relation to the matrix elements A^™' n _ 1; A^fl\, A m n £ l } 1 and A^m.n^ of A^ m,n "> entering in eq. (47). 
This results in 

Cra-ln-lCmn Cm-lnCmn-1 = dct(A^ ' ^ )C m _i n _i. mrl . (53) 

Using the explicit form of the derivatives in the determinant of (46) the cofactors can be expressed as derivatives 

C mn =det^ m - 1 -"- 1 ), 
C m - ln = -S x d C tA (m - 1 - n - 1 \ 
C m -m-i = 5J y det ^C™- 1 .™- 1 ) , (54) 
Cmn-! = Sy det A^" 1 '"- 1 ), 
C7 — Het 4( m ~ 2 >™- 2 ) 

This results in the recursion relation 

[My det A^" 1 '"- 1 )] det A ( - m ~ 1 ' n ~ 1 ^ - [S x det A^" 1 '"" 1 )] [<f„ det A^" 1 '"" 1 )] = det ^( m >™) det A^" 2 '"" 2 ) . (55) 
Raising the indices of the recursion relation by one this relation can be rewritten as 

5 X S V log detM (m '™ } ) = 5— . (56) 

The Toda lattice equation for the generating functional of the two-point function of the chUE follows immediately by 
applying this identity to (47). It is given by 

SJ y log Z m > n {x,y) = Anmx y — . (57) 

[Z m ; n (x,y)\ 2 



3. The Replica Limit of the Toda Lattice Equation 
The microscopic disconnected scalar susceptibility for the chUE is given by 

^^^Jfco^^ 108 ^^ (58) 
Using the Toda lattice equation above in this definition we find 
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where, explicitly, 



and [9] 



X {v) {x,y) = 4xyZ[';}(x,y)Z^l_ 1 (x,y), (59) 



Z[%,y) = —L-( X I V+I (x)l v (y) - yl v+1 {y)l v {x)), (60) 

y x 



Z^l_^x,y) = - ir — 1[ {xK v+1 {x)K v {y) - yK v+l (y)K u {x)). (61) 

y x 



The two point spectral correlation function is given by the discontinuity of the disconnected scalar susceptibility across 
the imaginary axis. The result is in exact agreement with the analytical result for the two-point function found in 
[57]. Again, the replica limit of the Toda lattice equation explains that the spectral correlation function comes from 
the product of a fermionic and a bosonic partition function. 

III. TODA LATTICE EQUATION FOR NON-HERMITIAN THEORIES 

In this section, we derive, in the limit of weak non-Hcrmiticity, the spectral density for non-Hermitian random 
matrix theories from the replica limit of the Toda lattice equation. We will consider the weak non-Hermiticity limit 
of two systems: 1) A disordered system in an imaginary vector potential (the Hatano-Nelson model [42]) which is 
in the universality class of the unitary ensemble. We will derive the known spectral density of this system [45,65] 
from the replica limit of the Toda lattice equation. 2) QCD at non zero baryon chemical potential which belongs 
to the universality class of the chiral unitary ensemble. We will derive the quenched microscopic spectral density 
of the non-Hermitian Dirac operator. This is a new result that we have checked against a high statistics numerical 
simulation. For a general discussion of non-Hermitian random matrix theories, we refer to a recent review by Fyodorov 
and Sommers [66]. 

A. Disordered System at Nonzero Imaginary Vector Potential 

The spectral density of the Hamiltonian for a disordered system in an imaginary vector potential was derived in 
[45,65] using the supersymmetric method. In [45] this result was obtained starting from an ensemble of complex 
random matrices, whereas in [65] the starting point was a nonlinear cr-model which is also applicable beyond the 
ergodic domain. This nonlinear a model is applicable to the Hatano-Nelson model [42] and belongs to the universality 
class of the (almost Hermitian) UE. In this section we analyze its partition function for n fermionic flavors. We find 
that it satisfies a Toda lattice equation enabling us to derive the known result [45,65] for the spectral density by 
means of the replica trick. The replicated bosonic partition function that enters in the replica limit of the Toda lattice 
equation is evaluated using the Ingham-Sicgcl integral [25]. 

1. The fermionic partition function 
The static limit of the partition function for n fermionic flavors is given by [65] (see also the discussion in [67]) 

*•<«•> - ^mw } /d*-***"**™". (62) 

where the integral is over U(2n)/[U(n) x U(n)]. The matrix S 3 is defined by 

J- ,) m 

and 1 is the n x n unit matrix. The parameter o determines the strength of the imaginary vector potential, while y 
is the mass of the fermionic flavors (i.e. the imaginary part of the argument of the resolvent). In order to evaluate 
the integral we choose the parameterization [4,6] 
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Q = TtS 3 T, (64) 

with 

T = V*( C ° S ° Shl9 a )v, (65) 
1 — sin 6 cos 1 



V = rj " . (66) 



and 

M ( 

Here, «i e J7(n) and u 2 G U(n)/U n (l) and = diag(#i, • • • , 9 n ). The integration measure is given by the product of 
the matrix elements of the off-diagonal blocks of 

ST = VdTT^ 1 V^ . (67) 

The calculation of the Jacobian of the transformation from the variables {8T\2, <5T2i} to the variables {<iFV^, <5cos#} 
is elementary and is given by 

J = ]j2cos6> fe ]J(cos 2 (9 fe -cos 2 6l ; ) 2 . (68) 

k k<l 

The angular integrations factorize from the partition function and result in the volume of the integration manifold 

K=vol 2 ([/(n))/vor([/(l)). (69) 

In terms of the variables 

A fc =cos 2 (9 fe -sin 2 (9 fc , (70) 

the partition function can thus be written as 

V n _ 

'o\(TT(n\\ I , ■ 

k k<l 



Z n {y;a) 



2«(n-i)-™ n !vol(?7(n)) 
This partition function can be written as a determinant of derivatives 



f ndA fc n^- A o 2 ^ 2Q2Ar ^ (A ^ i)+2yAr ^ Afc - (7i) 



where 



Z ^ a ) = 2 2"("-DiV"("-i)vol(E/(n)) det [dg +g 3i(y;a)] P , g =o,...,n-i (72) 



Z 1 {y- a) = 2 J' dXe - 2a2N ^ 2 -^+ 2 y N \ (73) 



The partition function (62) therefore satisfies the Toda lattice equation (cf. the discussion in section II C 2) 

1 fl2 , -7 I \ 8nZ„ +1 (y;a)Z„_ 1 (y;a) 

— <9 2 logZ„ (y;a) = — . 

N 2 y ir Z 2 {y;a) 



(74) 



2. The Spectral Density 

The spectral density is given by (z = x + iy) 



p(x, y- a) = Inn 1 1 A A. log Zn(z , z *, a) . ( 75 ) 



Taking the replica limit of the Toda lattice equation we thus find 
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where (cf. (73)) 



and, as we will show below, 



8N 2 

p(x,y;a) = — r Z 1 (y;a)Z- 1 (y;a) (76) 



Z^y- a) = 4e 2a2jv / cosh(2yNt) exp (-2a 2 Nt 2 ) dt (77) 
Jo 



Z. x {y- a) = ^i7^ CX p(-a 2 7V) exp (-^ ) ^ 



The spectral density is x- independent. This implies that if we change the non-Hermiticity parameter a there is only 
spectral flow perpendicular to the j/-axis. Therefore, the normalization C_i of Z_i has to be chosen such that the 
integral 

% r dy Zl (y-a)Z^(y-a) = 3 -^l e ^ (79) 

is independent of a. Normalizing the integrated spectral density to N per unit length we find 

CLi^e""' 2 . (80) 
Combining the previous five equations, we find the spectral density 

p(x, y; a) = ^0- exp I cosh(2yNt) exp (-2a 2 Nt 2 ) dt, (81) 

V2na V la J Jo 

in complete agreement with earlier work [45,65]. Not only do we obtain the exact analytical result by means of the 
replica limit, but we have also shown why the spectral density factorizes into a bosonic and a fermionic partition 
function. 

Before we derive Z_\ let us comment on the a-dependence of the normalization integral (79). The partition function 
Z\ was derived from the effective partition function (62). The commutator term arises because of the requirement 
of local gauge invariance [43,44,65,68] of the effective Lagrangian with the kinetic term included. The non-Hcrmitian 
random matrix theory in [45] does not have this property. Starting from this random matrix model [45], which will 
be introduced in the next section, the exponential prefactor in (77) would have been exp(iVa 2 ) instead of cxp(2./Va 2 ) 
and the normalization integral would have been automatically a-independent. 



3. Calculation of Z-\ using the Ingham- Siegel integral 

In this subsection we calculate Z_i for a non-Hcrmitian random matrix model using the Ingham-Siegel integral 
as proposed in [25]. The saddle point manifold is obtained automatically as a result of this integral. This should 
be contrasted with the standard approach which uses the Hubbard-Stratonovitch transformation and where it is 
necessary to deform the integration contour in order to be able to interchange integrations [69] . 

The partition function Z_i is defined by 

Z-^^^.JT^i^i^ + H-A Z + H ( +A )), (82) 
where the average (• • •) is over the probability distribution 

P{H,A)=e-^ H ' H+ ^ A ' A 1 (83) 

and the N x N matrices H and A are Hermitian, = H, and anti-Hermitian, = —A, respectively. The 
normalization constant will be chosen such that the limit e — > is finite. In order to write the inverse determinant as 
a convergent bosonic integral 
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det" 1 ^ * € H ^ z + H + A ^ = J d(j) ie t[e ^ k ^ +e ^ k ^ + ^ k(zSkl+Hkl+A ^ (84) 
the imaginary part of e has to be positive. We collect the bosonic fields into the positive definite 2x2 matrix 

N 

Qij = ( 85 ) 

fe=l 

If we introduce a mass matrix by (<7j are the Pauli matrices) 

C = xa x - ya 2 (86) 

the partition function can be written as 

Z- 1 {z,z*-a)= lim C t ( d<t> l € lTr ^ T+t ^-^ TlQa ^-^ TrQ ^ Q ' y \ (87) 

e— >0,iV— >oo J 

The four-boson terms can be linearized by means of a 5 function, [46,70] 

Here both F and Q are Hermitian. This results in 

Z_ 1 (z,z*;a)= lim C £ / ^dQdFe^(- f + e ^+^^+ lTr ( cT ^-^v TrQfflC?CTl -^ TrQff2C?CT2 . (89) 
Under the assumption that the integrals can be interchanged we obtain 

Z_ 1 (z,z*;a) = lim C e [ dQdFdet~ N (F - e y^FQ+^(fQ)-^Q^Q^-^Q°*Q°*. 

e^0,N^oo J 

The integral over F is an Ingham-Siegel integral of the second kind. For Im(e) < it is given by [25] 

J dFdct- N (e + F)e lTlQF = C N , p 8(Q)det N - p (Q) e - ieTlQ , 



(90) 



(91) 



where the integral is over pxp Hermitian matrices, Cm. p is an irrelevant constant, and 9(Q) denotes that Q is positive 
definite. After rescaling Q — ► NQ we find (N > 2) 

Z_ 1 (z,z*;a)= lim C e /dg0(Q)det Ar - 2 O e jJVTr (^ 1 +^ 2+e ) Q )-* Tr< 3 CTl ^ 1 -^ TrC ? CT2( 3 CT2 . (92) 

We calculate this integral in the weak non-Hermiticity limit and the microscopic limit, i.e. N — > oo with a 2 N and 
zA^ fixed. In this limit, the saddle point equation, given by 

{<?iQ? = 1, (93) 

fixes two of the four degrees of freedom in Q. In order to complete the evaluation of -Z-i, we will choose an explicit 
parametrization of Q and perform the integration over the two degrees of freedom not fixed by the saddle point 
equation exactly. Positive definite 2x2 Hermitian matrices can be parameterized as 

n _ t f e " cosn s ie 1 ^ sinh s \ , . 

' \ —ie~ 1 ^ sinh s e _u coshs J' ^ ' 

The Jacobian for the transformation from the Q variables to the variables of this parameterization is given by 

J = 4e 4t cosh s sinh s. (95) 

The saddle point condition (93) is given by 

t = 0, = 0, (96) 
13 



whereas the variable s and u are not determined by the saddle point equations. They parameterize the saddle point 
manifold. The Gaussian fluctuations in the t and <p variables about this saddle point manifold give rise to a factor 
l/(2sinhs). As final result we thus obtain 



/'OO r-OO 
du / ds COsh se a2N(1 - 2 C ° Sh2 s)-i2Ny a inhs+2iN e cosh „ cosh, 



(97) 



The integral over u results in a contribution ~ — loge. This leading order singularity in e will be canceled by the 
normalization constant C e and the remaining prefactors, 2C e , will be denoted by C-\. The integral over sinhs can 
then be performed by completing squares. This results in the Gaussian y dependence (78) 

c -i 2 A n„„„/ N v 2 



*-i(y;a) = — ^-exp(-a^)exp j . (98) 

The Toda lattice equation suggests that the partition function with n bosonic flavors can be expressed as a deter- 
minant of derivatives. This leads us to conjecture 



Z_ n (y;a) = C_„ det^+^.^y; a)] M= o,...,n-i (99) 

nNy 2 1 

- (100) 



a" 

where Z_i is given by (78). This conjecture will be discussed in a separate publication [52]. 



B. The Phase Quenched QCD Partition Function 



In this section we derive the microscopic spectral density of the quenched Euclidean QCD Dirac operator at nonzero 
baryon chemical potential. The calculation is performed in the weak non-Hermiticity limit, where only terms to up 
to second order in the chemical potential contribute to the effective action. In this subsection we follow the usual 
convention that the massless Dirac operator at zero chemical potential is anti-Hcrmitian. Therefore "weak non- 
Hcrmiticity" is actually "weak non-anti-Hermiticity" , but we will refrain from using this name. 

The generating function for the quenched spectral density at nonzero chemical potential is not the QCD partition 
function with n flavors, but rather the QCD partition function with n flavors and n conjugate flavors [71,72]. Since 
the fermion determinants of the flavors and the conjugate flavors are each others complex conjugate this partition 
function is also known as the phase quenched QCD partition function. For n flavors with mass z and n flavors with 
mass z* it is given by (in the sector of zero topological charge) 

Z n{ z, = (det™ ( ./ + / d +/ ) det" ( _J\ ^ ~ id z t " ) , z = X + iy, (101) 

where id is the covariant derivative and the average over gauge field configurations is weighted by the Yang-Mills 
action. For positive n the low-energy limit of this theory in the chirally broken phase is [73] 

Z(ZZ*-U) = - f dUe -^f^'TAU,B}[U\B\ + ^VTr(MU+MU^) mq 2 \ 

with' 1 ' 

B = S 3 and M=( Z Q ^Y (103) 



^^We stress that the combination (MU + MU^) enters the action and not the usual combination (M T {7 + MU^). The reason is 
that the phase quenched QCD action (101) has n flavors of mass z and other n flavors has mass z* . The mass term that follows 
from chiral invariance is given by Tx(MrlU + Mlr.U^) with Mrl the right-left handed mass matrix and Mlr the left-right 
handed mass matrix. The usual mass term is obtained if Mrl = M LR . We however have that Mrl = Mlr = M with M 
complex. 
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The volume of space-time is denoted by V. The resolvent is defined by 

G(z,z";fi) = lim \d z log Z n (z, z*; //). (104) 

At nonzero \x the support of the spectrum of the Dirac operator is a two-dimensional domain in the complex plane 
with spectral density given by 

p(z; t i) = -d z .G(z,z*;n). (105) 

In the thermodynamic limit at fixed fi and m, the effective partition function (102) can be analyzed in terms of mean 
field theory [68,73,74]. The result is [73] that the spectral density is constant inside the domain 

2u 2 F 2 

\Rez\ < -L (106) 

and zero outside this domain. 



1. Toda Lattice Equation 

Here we will analyze the effective partition function in the microscopic limit where both /j, 2 V and zV remain fixed 
in the thermodynamic limit. Our aim is to calculate the spectral density of the quenched theory. To do this we will 
show that Z n satisfies the Toda lattice equation. The proof follows directly from the proof for two point function at 
zero chemical potential: In the parameterization discussed in section II C the trace of the first term in the action is 
given by 



Tr[C/,B][C/t, jB ]=8^(A 2 fe -l). (107) 
fe=i 

Repeating the steps (43)- (46) with this factor included we find (we recall the notation S z = zd z ) 

Dn 

( zz *}n(n-l) 



Zn(z,z*;») = det [SfoZ&Xw)]^,^, (108) 



with 



Zi{z,z*\n) = -e 2VF ^ 2 f 1 d\\e- 2VF ^ 2x2 1 (\zZV)Io(\z*XV) 

K JO 



(109) 



and 

2n(n+l)/2 
7r n(n-l)/2 n «-l k l 

This results in the Toda lattice equation 



D n = —zr^^rr-,- (no) 



x x i r, , * n nn , ^2 Z n+i{z,z*;y)Z n ^ l {z,z*;y) 

6 Z S Z , log Z n (z,z ;n) = —{zz ) ^ -—^ . (Ill) 



The quenched spectral density is then given by 



p{x,y;n) = ]im —d z d z . log Z n (z,z*;n) 

= Hm zz* Z n+1 (z,z*;iJ,)Z n - 1 (z,z*;n) 
n™ 2 Z*(z,z*;n) ' 1 ' 

In our normalization Z (z, z*) — 1 so that 

zz* 

P{x,y,fj) = —Z 1 (z,z*;fi)Z- 1 (z,z*;n). (113) 

Our next task is to calculate Z_\{z, z*; fi). To do this we derive the static low energy partition function from a chiral 
random matrix model corresponding to (101). The main problem is to choose the correct integration manifold. As is 
discussed in the next section, this can be done without guess work if we use the Ingham-Siegel integral as proposed 
in [25]. As an introduction to this method we have included the calculation of Z_\(z, z*) for fi = in Appendix C. 



15 



2. The effective partition function Z-i 



In this section we derive the effective phase quenched partition function in the microscopic and weak non-Hermiticity 
limit. The integration manifold is determined by symmetries and convergence requirements. Usually, these conditions 
define a suitable integration contour for the cr-field. A method that avoids this somewhat ad hoc procedure was 
advocated in [46,25]. Below we will follow this approach in the calculation of the phase quenched partition function 
for one bosonic flavor. 

The partition function for one bosonic flavor and one conjugate bosonic flavor is given by the integral 



lim [ 

JV->oo J 



dWe 



iW> + fi z 



(Mr 



z* -iW + n 
~iW^ + n z* 



(114) 



If z is inside the domain of eigenvalues, the partition function potentially diverges. Therefore, we have to regularize 
the determinants. A suitable rcgularization is one where the determinant can be written as a convergent bosonic 
integral. This procedure, which is also known as hermitization [75], amounts to rewriting the product of the two 
determinants as 



dct~ 



lim det 

£^0 



z iW + /i 
iW^ + [i z 



(Mr 



e 
e 

z* -iW + n 
\ -iW* + fj, z* 



z* -iW + n 
■iWi + n z* 

z iW + ^\ 
iW^ +ii z 
e 
e 



(115) 



After rearranging the rows and columns inside the determinant, the partition function (114) can be written as 



Z-i = lim C, 



/ 



dWe 



I' 



b k d<pl exp[i 




/ 


e 


z 


iW + ii \ 






z* 


e 


iW-fi 








-n 


e -z* 


03 


V 





-iWt 


-fi -z e J 


\4>A 



(116) 



After averaging over W the limit e — > is potentially singular. Such singularity will be absorbed as a multiplicative 
renormalization constant into C e . As in the case of the Hatano-Nelson model the overall prefactor of Z_\ will be fixed 
by the normalization of the complex density. For this reason we do not keep explicitly track of the prefactor, it will 
be simply labeled by C e . 

The Gaussian integration over W results in a 4-boson interaction term given by 



exp 



-^TrQ!Q 2 



(117) 



where 



Qi 



<?2 



i>4 '03 01 ' 04 



(118) 



and we have used the notation 



<t>i 



N/2 



(119) 



Instead of the usual Hubbard-Stratonovitch transformation, we linearize the 4-boson interaction term with the help 
of matrix 5 functions. For Hermitian matrices Qi and Qi the 5- function can be represented as 



S(Qi - Qi) 



(27T) 



tt) 4 / 



dFe -iTrF(Q i -Q i); 



(120) 



where the integral is over Hermitian matrices F. This results in 



16 



exp 



J ' dQ x dQ 2 J ' d FdGe^- iF ^-^- iG ^-^-* Q ^. (121) 

The integral over the (j>k is uniformly convergent in F and G which allows us to interchange the order of these integrals. 
We finally obtain the partition function 

Z_ x = lim C £ [ dQ,dQ 2 [ dFdGe^- iF ^- iG ^-%^det-% ( C + T ^ nT ) , (122) 

e— >0,iV — >oo J J V ' i<T3 _ ^ / 



where we have used a block notation and 



The anti-symmetric matrix / is defined by 



c = ( ;* \ ) . (123) 



1 

-1 



(124) 



We further simplify this integral by changing integration variables according to F — > F — ( T and G — > G + I( T I and 
Q 4 -> NQ l /2 7 i = l,2. This results in 



Z_i = lim C, 

e->0,JV^OO 



J dQ x dQ 2 J dFdGe^-^^^+^^-^-^^dct-^ ( 6 +^ ^ ) , (125) 

where, for reasons of convergence we also left infinitesimal increments inside the determinant. Contrary to the 
example in Appendix C, we were not able to perform the F and G integrations analytically. However, in the weak 
non-Hcrmiticity limit, where /j, 2 N remains fixed in the limit N — ► oo, the determinant is given by 



(126) 



In the limit N — > oo, at fixed /U 2 ./V, the F and G variables in the n 2 N term can be replaced with the saddle point 
values of F and G at = 0. They are given by 

iQu —^=*Q2- (127) 



e + F ^" e + G 

The remaining integrals over F and G are Ingham-Siegel integrals (see (91)). They can be performed exactly resulting 
in 

Z_!(z,z*; M )= lim C e /rfQidQ 2 e(Qi)0(Q 2 )det^^ 2 (QiQ 2 )e Tr ^ cT ( Ql - /Q2/ '-* QlC?2 -^ 2Qlff3Q2CT3 ], (128) 

where we remind the reader that £ contains the regulator mass e. (Recall that is a matrix step function which is 
equal to unity if Q is Hermitian and positive definite and equal to zero otherwise). In the limit N — > oo, the integrals 
over the massive modes can be performed by a saddle point approximation. The saddle point equations are given by 

0r 1 - ( ?2 = 0, Q2 1 -Qi=0- (129) 

Both equations can be rewritten as 

Qi = Q2 1 , (130) 

and therefore only four of the modes, which we choose to be Q2, can be integrated out by a saddle point approximation. 
The quadratic fluctuations give rise to a factor 7r 2 /det Q\. The integral over the remaining modes has to be performed 
exactly. We thus arrive at the partition function 



= limCe / -^_e(Q 1 ) e ^[^C T Wi-'0r 1 0-*M 2 0i«T.0r 1 '3]. (131) 
^0 J det z Qi 

Comparing this result with the fcrmionic effective partition function (102) we can make the following identification 

N< — ► V, 

S < — > 1, (132) 

2 2 T7.2 

M < ► M F n ■ 

We will evaluate this integral explicitly in section IIIB4. Before doing so we now rederive this partition function 
based on the symmetries of the QCD partition function. 
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3. Symmetries of Z-\ 




In this section we show that the effective partition function (131) is completely determined by the symmetries of the 
underlying microscopic partition function. This implies that any theory with the same pattern of global symmetries 
and spontaneous symmetry breaking is described by the same low-energy effective theory. The only memory of the 
microscopic theory are the two constants, £ and F n , in (131). In particular, this means that the low-energy limit of 
QCD and the random matrix model discussed in the previous section are the same. The virtue of the random matrix 
approach is that the low-energy limit of the partition function can be derived directly from the microscopic theory. 

For C = and \i = the random matrix partition function (116) is invariant under Gl(2)/U(2) x U(2). Explicitly, 
the axial Gl(2)/U(2) invariance is realized by 

( <t> i \ 

with U\U A X = 1. (133) 

For reasons of convergence we can allow only symmetry transformations that do not alter the complex conjugation 
structure of the partition function. Therefore, 

U A = e H , with H f = H. (134) 

The vector U{2) symmetry is realized by 

with U V U V = 1. (135) 

Imposing these symmetries at nonzero mass and chemical potential requires that, in the block notation, 

Cl ft 1 ) (136) 




the mass and chemical potential are transformed as 



(2^U A ( 2 U A , ^2^U A ^ 2 U A \ (137) 



under axial transformations and as 



Ci — » UvCiUy 1 , (jbi — » Uv^iUy 1 , 



(2 - UvbUy 1 , M 2 - Uv^Uy 1 (138) 



under vector transformations. The matrix Qi introduced in (118) transforms as 

Qi - UaQiUa, Qi - UyQ^, (139) 

under axial and vector transformations, respectively. To first order in the mass matrix and second order in the 
chemical potential, we therefore can write down the following nontrivial invariants (Qi has the same transformation 
properties as (Qi) 

TrQiCf, TrQ^Cj, TrQ^Qf Vf- (140) 

Using that 

Ci = C, (2 = -1(1, Mi=M2 = Mc r 3, (141) 
we find the nontrivial invariant terms 

TrQiC T , -TtQ^I^I, ^TrQ^aQ^aa. (142) 
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The effective partition function should be invariant under the discrete symmetry, ( — > —1(1, of the microscopic 
partition function as well. Because the integration measure on positive definite Hermitian matrices, dQi/det 2 Qi, is 
invariant under Qi — ► Q^ 1 this is the case if the mass terms occur in the combination 

TrQxC T - TrQ^I( T I. (143) 

The numerical prefactors of the invariant terms are not fixed by the global symmetries of the random matrix model 
but follow from matching to the microscopic theory. 

At nonzero chemical potential the phase quenched QCD partition function with one pair of bosonic flavors has to be 
hermiticized as in the random matrix model (116) resulting in a Dirac operator with the structure given in (116) and 
a mass matrix as in (123). The partition function has the global axial and vector symmetries of (133) and (135) if the 
mass and chemical potential matrices are transformed according to (137) and (138) as well as the discrete symmetry 
( — > —1(1. The Goldstone manifold is parameterized by the axial tranformations and is therefore the coset of positive 
definite matrices [21,22] which transform according to (139). At zero chemical potential, the effective Lagrangian 
density to order p 2 respecting these symmetries is therefore given by 

L-i = tfTrdpQdnQ- 1 ~ ^(TrQC T ~ ^Q'^Cl)- (144) 

The phase of the mass term is determined by the condition that the integral over Q is convergent for e > 0. For 
real quark masses, z — z* = m, the mass term is minimized by the saddle-point solution Q — —i<j\. The vacuum 
energy is therefore given by 2to£V which identifies S as the chiral condensate. As is the case in the supersymmetric 
formulation of partially quenched chiral perturbation theory [14], the pion decay constant in L_i is taken to be the 
same as in the case of fermionic quarks. 

The chemical potential enters in the QCD partition function as an external vector field, 

V„ = fia 3 S v0 - (145) 

Therefore, the vector symmetry (135) can be promoted to a local symmetry by transforming V v as a nonabelian gauge 
field 

V v -» Uy X d v U v + Uy X V v U v . (146) 

The effective Lagrangian at nonzero chemical potential should respect this symmetry. This is achieved by replacing 
the derivative in the kinetic term by the covariant derivative 

d u Q -» V„Q = 8 V Q + [V V ,Q] (147) 

with V v given by (145). The low energy effective Lagrangian density is therefore given by 

L-i = ^TrV.QV.Q- 1 - i^(TrQ( T - TWT 1 /^/). (148) 

In a parameter domain where the fluctuations of the zero momentum modes are much larger than the fluctuations of 
the nonzero momentum modes the zero momentum part factorizes from the partition function [50] and is given by 

Z^(z,z*;^) - \imC e j ^-e(Q)e^ i ^^ Q - IQ ~ ll ^ F ^ Q ^ Q ~ 1 ^ , (149) 

where dQ/det 2 Q9(Q) is the integration measure on positive definite Hermitian matrices. Up to an overall constant, 
this result agrees with the effective random matrix partition function imposing the identification (132). 

4- Evaluation of Z _i 

We parameterize the positive definite matrix Qi in (131) as 

„ _ t ( e r coshs e ie sinhs \ , . 

ye _te siiihs e~ r coshs }' 

The Jacobian of this transformation is easily found to be 
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J = 4e 4t coshssinhs. (151) 
The range of the integration variables is given by 

re (—00,00), se (-00,00), te(— 00,00), #e(0,7r) . (152) 

Inserting this parameterization in (131) we find that the partition function is given by the following integral 

Z_i = lim C £ J drdsdtd<9 cosh s | sinhs | e 4 * (4Ksinhs ^ (153) 

The integral over r results in the modified Bessel Ko(2Ne cosh s cosh t) with leading singularity given by ~ — loge. 
This factor is absorbed in the normalization of the partition function. Then the integral over 9 gives a Bessel function. 
Introducing u = sinh s as new integration variable we obtain 

/oo />oo 
dt / duuJ (2Nu(x 2 cosh 2 t + y 2 sinh 2 t ) 1 / 2 )e~^ N ^ 1+2u2 \ (154) 
-oo JO 

where the finite prefactor has been labeled C_i. The integral over u is a known integral over a Bessel function. The 
final result for the phase quenched partition function with one bosonic flavor is thus given by 

C IP-^ f°° W^ 2 cosh 2 t + H 2 sinh 2 t) 

C _ ie -»» a m4^1 T , t N(x 2 + y 2 ) ^ 



5. The Spectral Density in the Complex Plane 
Reminding that, with the identification (132), the partition function for one fermionic flavor is given by 

Z 1 (z,z*;n) = -e 2N ^ dXXe~ 2N ^ x2 I (XzN)I {Xz*N), (156) 

I" JO 

the final result for the spectral density of the chGUE in the weak non-Hermitian limit obtained from the replica limit 
(112) of the Toda lattice equation (111) is given by [z = x + iy) 

ZZ 

P{x,y,n) = —Zi{z,z*;[i)Z- 1 (z,z*;n) 



8nfi 2 N 



■(x 2 + y 2 )e^^K ( N{x2 A + 2 y2) ) C d\Xe- 2N ^ x2 I a (XzN)I (Xz*N). (157) 

Jo 

As in section (III A) the constant C_i can be obtained from a spectral flow argument. In this case, the spectral flow 
with variation of fi becomes perpendicular to the y axis only for Ny 2 3> \i. In this limit we can use the large argument 
asymptotic approximation of the Bessel functions, 

Ko{z) = x IJ-e- z , / (z) = -i=( e z + Je - z ), h{z') = -;L={e*' -ie~*'), (158) 

resulting in the large y asymptotic limit of the spectral density (note that the oscillating terms in [z — z*)NX arc 
subleading after integration over A) 

y> A = 4^)3/^5/2/ ^" l dXXe-™^ 2 (e 2 ^ + e"™-). (159) 
The integral over x of this expression can be calculated analytically, 



j — i 



p(x, Ny 2 » ^ 2 ; rfdx = (160) 
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and should be independent of /j,. In the normalization of N eigenvalues per unit length in the y direction this results 
in 



C_i = 4wN 4 e- N » 2 . (161) 
The microscopic spectral density, obtained in the limit N — > oo with [i s = p?N and i] — zN fixed, is thus given by 

= ^ t^+I^ ^ r 1 dAAe -2^A 2 Jo(A )/o(A (162) 

2/^1 8pt^ / 

Again we note that Z_i was calculated starting from a chiral random matrix model, whereas Z\ was obtained from 
an effective partition function with a commutator term that was a remnant of the local gauge invariance of the chiral 
Lagrangian. If we also would have derived Z\ from the random matrix model, the exponential prefactor would have 
been e Nfl instead of e 2Nfl , and the normalization integral would have been automatically /i-independent. 
The leading order asymptotic expansion of the bosonic partition function in /j 2 /N(x 2 + y 2 ) is given by 

U " V ' W 2^N^2N(x 2 +y 2 ) K \N{x 2 +y 2 )) > V ; 

If we insert this result in the expression for the spectral density (157) we reproduce the result in [48] which was 
obtained by means of the complex orthogonal method [76] for a closely related partition function defined in an 
eigenvalue representation. The discrepancy between the phase quenched partition function (102) and the partition 
function in [48] is even more remarkable since it was shown [77] that, in versions of both models with the phase of 
the fermion determinant not quenched, the two models coincide in the weak non-Hcrmiticity limit. 



Exact Result 
Asymptotic Limit 



100 




— , | i i i r- 

Exact Result 
Asymptotic Limit 



0.02 



FIG. 1. The radial microscopic spectral density of the quenched QCD partition function at nonzero chemical potential. The 
histogram shows the result of a numerical simulation in which 2.000.000 random matrices of size 200 x 200 where diagonalized. 
The value of the chemical potential is given by /j, 2 N — 0.04 with N = 200. The full line is the exact analytical result and is in 
agreement with the data. Also shown is the result with the asymptotic limit of the bosonic partition function in the expression 
for the density (dotted curve). This curve fails to describe the data when r//tt is not in the asymptotic domain r/fi S> 1. The 
left figure is a blown-up version of the right figure for r < 0.005. 
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To convince the reader of the correctness of our result, (162), we have performed a numerical diagonalization of 
2000000 matrices of the form 



iW + p 
iW^ + p 



(164) 



The matrices where of size N = 200 and drawn on the weight (15). In Fig. 1, we have plotted a histogram of the 
radial density, 

p(r;p) = / d(j)p(x = r cos <p, y — r sin (f>; p) (165) 
Jo 

for p 2 N = 0.04. In the same figure we plot our analytical result (solid curve) and the leading order asymptotic 
approximation of our result (dotted curve) which coincides with the result of [48]. Recently, Dirac spectra were 
calculated for quenched lattice QCD at nonzero chemical potential [78] . Although the lattice spectra were shown [78] 
to be in agreement with [48], the statistical accuracy is not sufficient to distinguish between the two analytical results. 
We are looking forward to simulations that will resolve this issue. 

In the limit y 2 N ^> p 2 and p 2 N ^> 1 the eigenvalues are distributed homogeneously in a strip parallel to the y-axis 
as has been shown by a mean field argument [73]. In this limit the Bessel functions in (162) can be approximated by 
their asymptotic expansion (see (159)) resulting in 

p s (v, V*;»s) = — ^= cxp(~ {V + V * )2 ) / dXe-^ie^+vn+e-Mn+vny (166) 
2p s V2ir op% J 



The saddle point of the A-integral is located at 



If this saddle point is outside the integration domain, the integral is exponentially suppressed. We thus find that in 
the limit p s 3> 1 and Im(?7) 3> p 2 the spectral density vanishes for 

If the saddle point is inside the integration domain, a saddle point approximation gives 



„<,,,.;„) = jL ,„ t !2±£1-$<I, 



in agreement with the mean field result (106) [73] provided that we make the identification (132). The total number 
of eigenvalues in a strip of unit width perpendicular to the y-axis is thus given by 

N 2 p s ( V ,r 1 *;p s )4p 2 =N, (170) 

in agreement with the choice of our normalization constant in (161). 

Finally, let us stress that we expect that it is possible to express Z_„ as a determinant of matrix with elements 
given by derivatives of Z-\. This leads us to conjecture the following form of Z~ n 

z -n(z> **;**)= (zz ^ ( m n-i) det **; /*)] w =o,i,...,»-i • (171) 

We hope to prove this conjecture in a future publication [52]. 

IV. CONCLUSIONS AND OUTLOOK 

We have shown that the replica limit of the Toda lattice equation is a powerful tool to derive the spectral correlation 
function of both Hcrmitian and non-Hcrmitian random matrix theories. We have obtained one-point and two-point 
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functions of Hcrmitian random matrix ensembles in the class A (Wigner-Dyson) and in class AIII (chiral) and one- 
point functions for non-Hcrmitian random matrix theories also both in class A and in class AIII. In the case of the 
chiral ensemble, the non-Hcrmiticity arose because of a chemical potential that was introduced in the random matrix 
model consistent with the covariance properties of the QCD partition function. Therefore, this result is particularly 
relevant for the analysis of Dirac spectra of quenched QCD at nonzero chemical potential. We emphasize that our 
analytical results for this case were not previously known. These results convincingly show that it is possible to derive 
nonperturbative results by means of the replica trick contradicting the lore that the method only can be trusted for 
the calculation of asymptotic series of analytical expressions. In addition, the replica limit of the Toda lattice equation 
explains the factorization of correlation functions in terms of a product of a fermionic and a bosonic partition function. 
Such a structure is present both in Hermitian and non-Hcrmitian random matrix theories. 

Although the results for the correlation functions in the other cases we have considered were known previously 
several new insights have emerged. In particular, we have conjectured that the fermionic, the bosonic and the super- 
symmetric partition functions are r-functions of a single integrable hierarchy which are related by means of the Toda 
lattice equations. In fact, this is the basis for the success of the replica trick in this approach. We have shown the 
validity of this conjecture for the fermionic generating function of the chGUE two-point function and the fermionic 
generating functions of the non-Hermitian theories. In several other cases this property was known previously. In the 
bosonic case, only the generating function for the GUE two-point function was shown to satisfy the fermionic Toda 
lattice equation for an arbitrary number of replicas. Of course, in cases were the spectral correlation functions were 
known previously it is clear that the Toda lattice equation can be continued at least to the partition function with 
one bosonic flavor. The difficulty with bosonic partition functions is that the weight in the integral has poles due 
to the inverse determinant. For Hermitian random matrix theories the poles are on the real axis and can be easily 
avoided. In the case of non-Hermitian random matrix theories with eigenvalues scattered in the complex plane, the 
poles are located in the very domain where we wish to calculate the spectral density. These problems can be dealt 
with by hermiticizing the generating function and using the Ingham-Siegel integral to avoid ambiguities in the choice 
of the integration contours. Although we have only evaluated the bosonic partition function for one bosonic flavor, 
we do not expect major technical problems for an arbitrary number of flavors and we expect to find the result that 
we have conjectured based on the Toda lattice equation. This problem will be addressed in a future publication. 

Recently, quenched lattice QCD Dirac spectra at nonzero chemical potential were compared to a non-Hcrmitian 
chiral random matrix model that was formulated in terms of eigenvalues so that the method of complex orthogonal 
polynomials could be applied to this model. The spectral density of this model differs from our result, but since 
the leading order asymptotic expansion of the two models is the same, the two results are close and, within the 
statistical accuracy, they are both in agreement with the existing lattice data. However, in a numerical simulation of 
a chiral random matrix model at nonzero chemical potential, we convincingly have shown that the model of Akemann 
disagrees with the data, whereas our result is right on the mark. Based on universality arguments one would expect 
that the two models would coincide in the weak non-Hcrmiticity limit, but apparently, the model of Akemann is in a 
different universality class. 

Having shown that it is possible to obtain exact analytical results by means of the replica trick some comments 
on the critique of the replica trick set forward in [4,7] are in order. One of the points that were raised is that the 
fermionic partition function or the bosonic partition functions alone are not sufficient to reproduce an exact analytical 
result. This is particularly clear in the fermionic case where the resolvent is an analytic function of the mass for all 
positive integer values of the replica index, so that its discontinuity is zero even in the replica limit. A careful analysis 
[7] shows that one cannot expect to obtain the full analytical answer in the bosonic case either. The conclusion of [7] 
was that, in order to derive nonperturbative results, the proper procedure is to combine information from bosonic and 
fermionic replicas. This is exactly what is achieved by the replica limit of the Toda lattice equation. However, the 
Toda lattice equation itself does not justify an analytical continuation from integer values of the number of replicated 
flavors. Only in the case of the one point function of the unitary ensemble with finite size matrices is it known how 
to derive a Toda lattice equation for any real number of flavors. 

It would be of great conceptual and practical interest to generalize the results obtained in this paper to the 
orthogonal and symplectic ensembles. The orthogonal chiral ensemble is relevant for QCD with two colors where 
lattice simulations at nonzero chemical potential already have been performed. In this case there is no sign problem 
and one can go beyond the quenched approximation. Several predictions for macroscopic properties [44,68,73,79] 
have been confirmed by means of lattice gauge simulations [80] . A comparison to the microscopic spectral correlation 
functions will offer an independent crosscheck of both numerical and analytical methods. 
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Appendix A. NOTATION 

In this appendix we comment on the notation used in this paper. In the main text we consider several the generating 
functions or partition functions. A partition function in the class of the Unitary Ensemble (abbreviated UE) will be 
denoted by Z and a partition function in the class of the chiral Unitary Ensemble (abbreviated chUE) by Z . If the 
probability distribution of the matrix elements of matrices in these ensembles is Gaussian, they will be called the 
Gaussian Unitary Ensemble (GUE) and the chiral Gaussian Unitary Ensemble (chGUE), respectively. If the form of 
the partition function does not depend on the specific form of the probability distribution we often omit the Gaussian 
adjective. The subscripts are used in the same way for both ensembles. We denote the partition function for m 
flavors of mass x and n flavors of mass y by Z m>n {x,y) and Zm]n{x,y). The additional superscript v of the chUE 
partition function denotes the topological sector. This superscript will be dropped for the trivial topological sector 
(i.e., v = 0). The sign of negative subscripts will be written out explicitly. It indicates that the corresponding flavors 
obey bosonic statistics. At nonzero imaginary vector potential or chemical potential the partition function depends 
both on the complex mass z and it complex conjugate mass z*. We use the notation z = x + iy and thus have 
Z n (z, z*;a) — Z n (x, y; a) and Z%(z, z*; — Z^(x, y; /i), respectively. Following the usual convention in the literature, 
the operators or matrices will be chosen Hermitian or weak non-Hermitian in the GUE class and anti-Hcrmitian or 
weak non-anti-Hermitian in the chUE class. However, for both case we will use the term "weak non-Hermiticity" 
when appropriate. 

Throughout this paper we will use the unnormalized Haar measure. 



Appendix B. CALCULATION OF A JACOBIAN 

In this Appendix, we calculate the Jacobian of the transformation 

mi \ fvi \ . / v\ 



u =v »a «rr 4)- ,172) 

where U is an (n + m) x (n + to) unitary matrix and A is the block diagonal matrix given by 



Ki={^ 1 ^ ft 2 1 > M = l,"-,2ro, 

A k j = 0, k > 2to or I > 2m, k ^ I, 

A Kk = -1, k> 2m. (173) 



The diagonal matrix fj, — diag(/zi, • • • , fj, m ), an d we will use the notation that = y/l — ^| with Afc € [0, 1]. Here, u\ 
and v\ are unitary to x to matrices, U2 G U(n) and V2 € U (n) / (U m (1) x U(n — to)). One easily verifies that the total 
number of parameters on both sides of (172) is the same. For the calculation of the Jacobian of this transformation 
it is convenient to use the differentials 



su 


= ( 4 4> 


8u\ 


= v\u 1 1 du\Vi, 


8u 2 


A -I 

= V^U 2 dU2V2, 


Svi 


— t>i dv\, 


Sv 2 


= v 2 ~ 1 dv2 


SX 


= dX. 



U^dU 



Vl 

V2 



(174) 

The matrices SU, 8uk and 5v k are anti-Hermitian. To calculate the Jacobian we have to distinguish six different cases 
depending on whether SU is diagonal or cither of its indices are smaller or larger than 2m. We split SU into 4 blocks 
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5U = 



SU U SU 12 
5U21 61/22 



(175) 



of size to x to, to x n, n x m and nx n, respectively. 

The Jacobian of the transformation from {SU^l, SU12, SUlfl, 5U 22 } to the variables {5v\ l , 6v 2 l , Su kl , Su* 1 , SXk} will be 
denoted by Ju- For 1 < k < to and 1 < I < to and k ^ / the Jacobian of transformation from {8U k {, SUi^SU^^SU^} 
to {5v± l , Sv* 1 , Su hl , Su^ 1 } is given by the determinant 



Jk< 



m,l<m,k^l 



det 



/ X k X t - 1 



— /XfeA; Xk^l — Hk^i 
— \k^l MfeA; —XkUl 
XkXi — 1 iik^i XkXi 



xf-x 



k- 



(176) 



where we have used that /x| + X\ = 1. The Jacobian of the transformation from {SU^^SU^} to the variables 
{5v2 l , SU2 1 } for 1 < k < m and to < I < n is given by the determinant 



Jk<m,i> m = det 
= Mfc- 



[i-k H>k 

-Afe - 1 -Afe 



(177) 



The Jacobian of the transformation from {SU^l, SU^i} to the variables {Sv^ 1 , Su* 1 } for k > to and 1 < Z < to is given 
by the determinant 



>m,l<m 



det 



Mi Mi 
-A;-l -A, 



(178) 



The Jacobian of the transformation from {SU^j} to the variables {Su^ } for k > to and Z > m is simply given by the 
determinant 



Jk>r 



= 1. 



(179) 



Finally, we consider the Jacobian of the transformation of the diagonal matrix elements. For the transformation 
{SU kk , 511*2, 6U21, SU22} to the variables {Sv kk , Su kk , 6u kk , SX k } with k < to we find the Jacobian 



/ 



rn , k < m 



det 



2Xk 



Vk 



A 2 fe 







AfeMfc MfcAfe —Xk/J-k —1/lJ-k 
MfeAfc Afe^ife —fikXi 1/nk 




\ 



/4 



Al 



(180) 



For k > m the only nonzero derivatives are (5J7|^ / Suf, k = 1 resulting into a Jacobian of unity. 

Multiplying the Jacobians calculated in previous paragraph we find the Jacobian of the transformation (172): 



J= LI (A 2 fe -Af) 2 H(2A fc )M 



2(n— m) 



(181) 



l<k<Km 



k=l 



Appendix C. CALCULATION OF Z_i FOR ^ = 

This Appendix illustrates the use of the Ingham-Siegel integral for the chiral unitary ensemble at = which was 
already considered in [25]. In this case there is no need to hermiticize the Dirac operator, and the inverse determinants 
can be represented as 
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det_1 ' zT^t T) det_1 (-^t f 



/ / n" ^ ' 



Other constants that enter in the calculation will also be absorbed in the irrelevant normalization constant C. The 
integrals converge if Re(z) > 0. The Gaussian average with probability distribution as in (114) results in (up to an 
overall constant) 



exp 

with 



(183) 



Ql -\tt-<h. tt-fo)' 0a " «-*J' (184) 

Introducing (^-functions as in (120) and performing the ^-integrations we obtain 

Z_i = C J dQ lC lQ 2 J dFdGe^ iFQ i +iGQ *-< Q i-< Q *-* Q ^det- N '\F)det- N / 2 (G), (185) 
where an infinitesimal imaginary increment is included in F and G and £ is now defined by 

The integrals over F and G can be calculated by means of the Ingham-Siegel integral (91). After also rescaling 
Q t -> NQ t /2 we find 

z -i= C J dQ 1 dQ 2 6(Q 1 )9(Q 2 )e T ^- N ^- N ^-^^dct^ 2 (Q 1 Q 2 ). (187) 

In the microscopic limit where £N is fixed as N — > oo the saddle point equations are given by 

QiQ 2 = 1. (188) 

We choose the Q\ variables to parameterize the saddle point manifold. The Gaussian fluctuations about Q 2 = Q^ 1 
give rise to a factor l/det 2 (Qi). As final result we obtain 

Z_ 1=C [ 6»(Qi)e TrhAr ^ (Ql+Q i _1)1 . 
J det Qi 

We parameterize the positive Hermitian matrices as 



(189) 



Qi = U{ ^ e i ) U~\ (190) 



The integration measure is given by 

dQ 



det 2 Qi 



— = (e Sl - e S2 )(e- Sl - e- S2 )d Sl ds 2 dU. (191) 



The integral over U can be calculated as an Itzykson-Zuber integral [81]. This results in 



oo 



g—Nz cosh si — Nz*' cosh s-2 ^ — Nz cosh S2—Nz* cosh si 



Z -^= C I dS ^l 3 *»^-e»){e-«-e-«) ^ _ z * )(coshsi _ coshs2 ) ■ (««) 



Writing z = x + iy, the limit y — > of this partition function is given by 

d*i / ^(e 51 - e S2 )( e - sl - e-Mje-JVxCcoehai+codia:,) = 8 c(^2(iv x ) _ if 2 (Nx)), (193) 

-oo J —oo 

which, up to an (arbitrary) normalization constant, agrees with the result obtained by means of different methods 
[9]. For z purely imaginary we did not succeed to further simplify the integral (192). 
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Appendix D. THE fj, = LIMIT OF EQ. (131) 



To convince the reader that the partition function with one bosonic flavor is correct we consider in this Appendix 
the \i = limit of (131) which is given by 

Z_!(z,z*; M ) =C [ _^-fl(Q 1 ) e ^[^C T (0i-/0r 1 0]. (194) 
J det Qi 



where 



1 \ , , / 



1= I _! o and C= U* ]■ (195) 



As in the previous Appendix, we parametrize positive definite Hermitian matrices as 



Qi = C/( e ei)^ 1 (196) 

with integration measure given by 

^ = (e Sl - e S2 )(e- Sl - e- S2 )d Sl ds 2 dU. (197) 
dct 2 Qx 

To rewrite the integral over U as an Itzykson-Zuber integral, we use the property of 2 x 2 matrices 



so that 



TrC T (Qi - IQ^I) = TrC T Qi + ^q~J<Qi 

= K(f + ^-OQ, (199) 
Evaluating the Itzykson-Zuber integral we find as final result 

Z - 1 - C /..* 1 /..*»<«* - «*" e "" -'-"> i( iV/2)(A,-A 2 )(a..-e..) ' < 200 » 
where Ai and A2 are the solutions of 

A 2 = (z + z*e- Sl - S2 )(z* + ze- Sl - S2 ). (201) 

We did not succeed to show that the partition function (200) agrees with (192) for arbitrary complex values of z. 
The reason is that the integration contour has to be modified appropriately so that the integral becomes convergent. 
However, for z purely real or purely imaginary, (192) can be recovered form (200). For z = iy (with y real) one easily 
shows that the integrands of the two partition functions are the same. This is not the case for z — x (with x real). In 
this case a convergent integral is obtained if we first shift the integration contours in (200) by 

7T 7T 

si — ► si + -i, s 2 —> s 2 - —i (202) 



in the first term and by 



in the second term. Wc then find 



7T 7T 

•si -» si - -i, s 2 -> s 2 + —i (203) 



/OC f-OC 
dsi / ds 2 (e Sl +e S2 )(e- Sl + e^ 2 ) 
-00 J — 00 



g— A^o:(cosh si+cosh S2) _|_ g — ^^(cosh S2+cosh si ) 

— 2A^x(coshsi + coshs 2 ) 



(204) 
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and finally obtain 

d x [xZ-!(x)] = 8{K$(Nx) + Kl(Nx)). (205) 
The solution of this equation is given by 

Z_y(x) = 8C(K$(Nx) - Kl{Nx)), (206) 

in complete agreement with (193). 
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